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We discuss the roles of the macroscopic hmit and of different system-environment interactions in 
the quantum-classical transition for a chaotic system. We consider the kicked harmonic oscillator 
subject to reservoirs that correspond in the classical case to purely dissipative or purely diffusive 
behavior, in a situation that can be implemented in ion trap experiments. In the dissipative case, 
we derive an expression for the time at which quantum and classical predictions become different 
(breaking time) and show that a complete quantum-classical correspondence is not possible in the 
chaotic regime. For the diffusive environment we estimate the minimum value of the diffusion 
coefficient necessary to retrieve the classical limit and also show numerical evidence that, for diffusion 
below this threshold, the breaking time behaves, essentially, as in the case of the system without a 
reservoir. 
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I. INTRODUCTION 

The problem of understanding the classical world from 
quantum theory is subtle and especially challenging when 
dealing with classically chaotic systems. 

Even the definition of classical chaos cannot be directly 
translated to quantum mechanics. Indeed, exponential 
sensitivity to initial conditions, used to define classical 
chaos, relies on the concept of individual trajectories 
in phase space, which is absent in quantum formalism. 
The use of classical phase space distributions, instead 
of trajectories, seems to be the way to circumvent this 
problem, since they can be readily compared with quasi- 
probability distributions defined for the corresponding 
quantum system. 

One expects, however, that the dynamics of the quan- 
tum and the corresponding classical system should differ, 
after some time, even if the initial distributions coincide. 
This time, often called Ehrenfest time or breaking time, 
while large for integrable systems, can be very short for 
chaotic systems. Indeed, in this case it has been shown Ij 
to be proportional to the logarithm of the inverse of an 
effective Planck constant, HcS, which is the ratio between 
Planck constant and a typical action of the system. For 
integrable systems, on the other hand, it scales as an in- 
verse power of Soft. In fact, quantum corrections become 
important when the distribution is able to explore the 
non-linearities of the potential, which can occur in a log- 
arithmic time scale due to the exponentially fast stretch- 
ing of the distribution, imposed by chaotic dynamics. 

One way to face this problem is to go to the macro- 
scopic limit, namely heff — > 0, resulting in an infinite 
breaking time Tfi,. Nevertheless, for any physical system, 
fie// is not zero and therefore t/j has a finite value, which 
can be short, even for macroscopic systems. Indeed, it 
has been argued that, due to the shortness of the sepa- 
ration time, even components of the solar system would 



exhibit quantum features, which is in contradiction with 
observation 

Reconciliation of quantum and classical predictions 
in this case is provided by the irreversible coupling of 
the system with an environment 0, 0, Q, which leads 
to the elimination of the quantum signatures, so that 
quantum and classical evolutions remain alike. In sys- 
tems that, isolated, exhibit dynamical localization, it 
was shown (Tj, ||8j, |3, ilfll] that noise and dissipation can 
strongly alter the situation and under certain conditions 
restore the classical-like momentum diffusion. For partic- 
ular choices of environment, it has been also shown that 
this reconciliation is possible under some conditions that 
involve a scaling relation between the effective Planck 
constant, the non-linearity parameter and the strength of 
the system-reservoir interaction [lIlollTlL llllTaL ITlllS j . 

One of the aims of this paper is to further explore 
the conjoint role of the macroscopic limit (fie// ^ 1) 
and of the interaction with the environment in the 
quantum-classical correspondence. In particular, we are 
also interested in the regions where the conditions for 
classicality do not hold and examine the time scales 
at which the quantum-classical correspondence breaks 
down. We should emphasize that this time scale is dif- 
ferent from the one present in the localization problem 
studied in 7, 8, 9, lOj since this phenomenom can be ab- 
sent in the model we are considering Another goal 
of this paper is to analyze in detail the impact of different 
system-environment couplings within the framework of a 
model that can be implemented experimentally. 

For this purpose, we revisit the kicked harmonic oscilla- 
tor (KHO), which has been the subject of studies both in 
the classical 17] and the quantum description [T^ITMIl^ . 
Despite some peculiarities and numerical difficulties pre- 
sented by the KHO, the possibilities of implementation 
with current available technology in ion traps [20] turn 
this model into a very attractive one. Moreover, in ion 
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traps, one is able to create artificial reservoirs |2l| and, 
in particular, different kinds of system-environment inter- 
actions have already been produced experimentally [2^ . 
This favorable scenario becomes complete with the possi- 
bility of tuning the effective Planck constant by changing 
experimentally accessible parameters. Indeed, heff and 
its scaling properties are related to the so-called Lamb- 
Dicke parameter (to be defined later in this paper) , which 
can be modified either by changing the trap frequency or 
the directions of the laser beams interacting with the ion. 

We analyze two limiting cases of environment coupling: 
zero temperature, which leads in the classical limit to dis- 
sipation without diffusion, and a reservoir that leads to 
diffusion without dissipation. This is not an unrealis- 
tic situation: indeed, the first may be mimicked by the 
sideband cooling mechanism in ion traps, under proper 
conditions, while the second corresponds to a white-noise 
position-independent random force, coupled to the oscil- 
lator. This is known to be the most important source 
of decoherence in actual experiments [23|. We present 
analytical and numerical results concerning the "dis- 
tance" between quantum and classical predictions and 
the breaking time. More specifically, in the dissipative 
case, expressions for the breaking time in three different 
parameter regions are derived and their physical conse- 
quences are discussed. In particular, in the most inter- 
esting region, we have a result similar to the one ob- 
tained recently, using another method, by lomin and 
Zaslavsky 121^ . For the diffusive environment we estab- 
lish the minimum value of the diffusive constant in order 
to restore the classical predictions, and we provide nu- 
merical evidence that the breaking time behaves like the 
one for the system without reservoir if the diffusion con- 
stant is kept below this threshold. We also show that the 
purely diffusive reservoir has a much stronger impact on 
the quantum-classical correspondence than the dissipa- 
tive one. 

In Section ^ we present the main features of the clas- 
sical model, both in the absence and the presence of cou- 
pling with the environment. In Section [llll we introduce 
the quantum model, discussing its connections with the 
experimental realization in ion traps and with the corre- 
sponding classical model. Section lTvl is divided into three 
parts, showing the results for the system without a reser- 
voir or interacting with dissipative or diffusive reservoirs. 
Two appendices contain detailed derivations of some of 
the results presented in the main text. 



II. CLASSICAL DYNAMICS: THE 5-KICKED 
HARMONIC OSCILLATOR 

The classical 5-kicked harmonic oscillator has been 
studied for both the isolated (without reservoir) and 
dissipative (2^ cases. Here we review the basics features 
of these models and present also the effects of the inter- 
action with a diffusive environment. 



A. System without reservoir 

We consider a particle of mass to in a harmonic poten- 
tial submitted to a sequence of periodically applied (5-like 
pulses. The Hamiltonian that describes this situation is 

222 

H = ^ + ^^^^^^+Acosikx)YS{t-nT), (1) 
2to 2 ^-^ 

n 

where v is the oscillator frequency, r the interval between 
two consecutive kicks and A their amplitude. The kicking 
potential is position dependent with a periodicity given 
by the wave vector k. 

The differential equations of motion due to this Hamil- 
tonian can be replaced by a discrete map. Between two 
kicks the system evolves accordingly to 

i + v'^x = 0, (2) 

while at the kicking times rtr there is just a shift in the 
momentum, so that 

= Xn, Pn ^ Pn + Aksin{kx), (3) 

where the variables immediately after and before a kick 
are indicated, respectively, by the presence or absence of 
the -|- superscript. After this integration we can connect 
the solutions before each kick by the following map 

Xn+i = cos(i^T)a;„ + sn\{i'T) / m v [p„ + A fc sin(/ca::„)] , 

(4a) 

Pn+\ = —mv sm{h'T)xn + cos(:/t) [pn + Ak sm(kxn)] ■ 

(4b) 

Using dimensionless variables v and u defined by 

V = kx , 

u — kp/mv , (5) 

the map becomes 

Vn+i — cos(q;)u„ -I- sin(Q;) [u„ -f K shi{vn)] , (6a) 



Un+i = — sin(Q;)i;„ -|- cos(q;) [u„ -|- K sin(vn)] , (6b) 

where K — Ak'^/mi/ and a = i^t are, respectively, the 
renormalized kicking strength and the ratio between the 
period of the kicks and the period of the oscillator. The 
system's phase space is unbounded and mixed, exhibit- 
ing stable islands surrounded by a stochastic web along 
which the system diffuses. The web is characterized by 
its thickness that broadens (shrinks) as the value of the 
chaoticity parameter K increases (decreases). For a = 
2TT/q {q integer), the stochastic web displays a crystal 
(g e (7c = {3, 4, 6}) or quasi-crystal symmetry((7 Qc). 
These basic features can be seen in Fig. ^ where a stro- 
boscopic plot for the map © is shown for q — 6 and 
K = 2.0 for different initial conditions. 
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FIG. 1; Stroboscopic plot for q = & and AT = 2.0, showing 
some stable islands as well as the stochastic web forming a 
hexagonal symmetry in the unbounded phase space. 

Another useful tool to study the dynamics is to follow 
the time evolution of the phase-space probability distri- 
bution. This method is especially suitable to problems 
where the notion of a single deterministic trajectory is 
absent as in the case of noisy and quantum dynamics. 
For the numerical evaluation of the density dynamics, 
one usually evolves an ensemble of trajectories generated 
accordingly to the initial distribution and then, by count- 
ing the fraction of trajectories lying in each cell of phase 
space, one recovers the density at a given time. How- 
ever, this method presents some drawbacks in our case, 
due to the unboundedness of the phase space and the 
consequent escape of trajectories. One way to get rid of 
this problem is to extend the phase space boundaries to 
be sure that, for the time scale one wants to simulate, 
no trajectories are lost. Nevertheless, the increase of the 
phase space area, keeping the size of the cells constant, 
requires a larger number of trajectories in order to get 
good statistics. This point imposes severe constraints for 
an efficient numerical implementation. 

Alternatively, we start with a uniformly distributed 
ensemble of trajectories, each of them carrying its own 
weight related to the initial distribution. Accordingly to 
this, the probability at each point is obtained by requiring 
that 

Pn{Vn,Un) = Po(wo,1to), (7) 

which means that the value of the initial probability Pq a-t 
every point (uq, wq) in phase space is transported to the 
image (w„, m„) of this point under ^he action of map © 
after n iterations. All the classical quantities calculated 
throughout the paper are obtained through the evalua- 
tion of each individual trajectory and then averaging over 
all of them taking into account the respective probabili- 
ties. It is important to mention that, although differences 
between trajectory-based and truly density evolutions are 
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FIG. 2: (Color online) Classical probability distribution for 
g = 6 and if = 2.0 after (a), 3(b), 6(c) and 9(d) kicks. 
The same structure of the trajectory-based stroboscopic map 
(Fig. is shown. After 9 kicks only the central structure 
is visible. The whole web structure would be seen for larger 
number of kicks. 



expected |27j , our simulations show no difference be- 
tween the two methods for sufficiently small phase-space 
partitions. 

Fig. 12 displays the evolution of an initial Gaussian 
probability distribution centered at the origin for the 
same parameters of Fig. ^ The numerical procedure to 
plot the distributions is similar to the one described to 
calculate the averages but, in order not to have problems 
with the dispersion of trajectories coming from neighbor 
regions in phase space at i = 0, we do the calculation 
backwards in time, choosing the grid at any instant of 
time t and evolving the points using the inverse map to 
find the probability of the inverse image of this point at 
t = 0. 

From Fig.|21one can identify the inner structure already 
present in Fig. ^ The whole web structure would be 
visible for larger number of kicks. 



B. The dissipative case 

The dynamics of the kicked oscillator changes if a dis- 
sipation mechanism is introduced. In our model this can 
be achieved by modifying the equatiqn of motion between 
the kicks (0) with the addition of a friction term propor- 
tional to the velocity 

i + v'^x + Tx = Q, (8) 
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where F is the dissipation rate. The discrete map ob- 
tained by the integration of © and the use of the shift 
in momentum ^ is 



-rr/2 



cos(a)a;„ 



sin(a) 



[pn + Ak sm{kxn)] I , 
(9a) 



Pn+l = e 



^^^1 — mfl sin(a)a;„ + cos(a) 

X [pn + Ak sin(fca::„)] |, (9b) 



where 



n ^ ^j,2 _r2/4, 

Pn/lTl = X + Tx/2 , 

a = fir . 



(10) 



We perform again a change to dimensionless variables 
v' and u', so that now 



w' = kx , 

u' — kp/mfl , 



(11) 



and 

w^+i = e"^""/^ {cos(aX + sin(a) [< + if' sin(t;^)]} , 

(12a) 

{- sin(aX + cos(a) [< + sm(w:j]} , 

(12b) 

where K' = Ak^/mfl. 

There are new scenarios arising from the addition of 
dissipation depending on the values of K' and Ft. When 
one of these parameters is changed, the system may 
change from a periodic to a chaotic motion in a sequence 
of period-doubling bifurcations In Fig. [3|we show 

this sequence and also the average Lyapunov exponent for 
K' = 6.0 and Ft/2 varying from to 1. For the bifurca- 
tion diagram we iterated the map H12|) for 10^ steps and 
plotted the last lO'^ points corresponding to the u variable 
on the vertical axis. The Lyapunov exponent is averaged 
over lO** different initial conditions, equally distributed 
around the origin, taking into account the same initial 
probability distribution as in Fig. [21 For each trajectory 
the exponent is calculated using the procedure described 
in [2^ for 10^ iterations. One can clearly identify the 
periodic regions corresponding to non-positive Lyapunov 
exponents and the chaotic ones where the system goes to 
strange attractors as the one shown in Fig. ^ 



C. The diffusive case 

Another kind of external disturbance that can affect 
the oscillator dynamics is the diffusion generated by noise 
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FIG. 3: Average Lyapunov exponent A (solid line) and bifur- 
cation diagram (dots) as a function of the dissipation parame- 
ter rr/2 for K' = 6.0 and q — 6. The horizontal line at u — 
was plotted as a reference for the Lyapunov exponent. For 
the bifurcation diagram, the vertical axis corresponds to the 
the last 10^ points of u after 10® iterations of the map 1121 . 
We only show the region from —2 to 2 for better visualization. 
The arrow refers to the case plotted in Fig. |3] 
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FIG. 4: Strange attractor obtained from the map 11211 for q = 
6, K' — 6.0 and rr/2 — 0.36 corresponding to a Lyapunov 
exponent of 0.697 in Fig |2l (denoted by an arrow there). 



due, for example, to fluctuating forces acting on the sys- 
tem. One possible mathematical description for the evo- 
lution of the probability distribution between two con- 
secutive kicks is the Fokker-Planck equation 



dP 
'dt 



dP 
ov 



dP 
du 



V 



Q2p Q2p 



(13) 



The first two terms describe the harmonic evolution while 
the third accounts for diffusion, both in v and u, with dif- 
fusion coefficient V. This diffusive term has two different 
effects on the system's dynamics. First, the noise limits 
the development of small-scale structures in phase space 
generated by the nonlinear dynamics, smoothing out the 
probability distribution. While the stretching of the dis- 
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tribution tends to generate thin structures in phase space, 
noise wiU tend to counterbalance this effect, leading to 
a lower limit in the width of such structures, which will 
depend on both the nonlinearity parameter and the noise 
strength 12j . The second effect is a faster spread of the 
system over phase space. In fact, the diffusion produced 
by the noise adds a new mechanism for connecting differ- 
ent parts of the web, thus enhancing the original chaotic 
diffusion. 



III. QUANTUM DYNAMICS: QUANTUM 
5-KICKED HARMONIC OSCILLATOR 

A. System without reservoir 

The quantum Hamiltonian for the (5-kicked harmonic 
oscillator is given by replacing in Eq. the variables x 
and p by operators: 



H 



P 

2m 



OO 



hKq COS (kx) ^{^ ~ ^'^) : 



(14) 



where we have defined Kg = A/H. 

It was shown in (2(| that this model describes the 
center-of-mass motion of an ion in a one-dimensional har- 
monic trap submitted to a sequence of standing-wave 
laser pulses, off-resonance with respect to a transition 
between the ground state and another electronic level. 
In this off-resonance condition the excited state is negli- 
gibly populated and can be eliminated adiabatically. The 
result of this elimination is an equation describing just 
the motional dynamics subject to recoils associated to 
the incoming laser pulses. 

In terms of the annihilation and creation operators a 
and for the harmonic oscillator, 



varying the direction of the kicking laser pulses with re- 
spect to the trap axis, since k stands, in this case, for the 
projection of the lasers wave vectors on the trap axis. 
In terms of "q, we can write: 



(a) — —(v + iu)=(v + iu), 
2r] 



K 



(18) 
(19) 



where (O, ltT31) and ((T7|) were used. While the classicality 
parameter 77 appears naturally in the quantum model, 
when this is expressed in terms of the annihilation and 
creation operators, it may also be introduced classically, 
by using the new variables v and u, which yields the 
following scaled map 



Vn+i — cos(a)w„ + sin(a) 



Un + TT sm(2?7w„) 
2r/ 



, (20a) 



Un+l 



' sin(Q;)w„ -I- cos(a) 



K 



Wn + sin(2?7u„) 
2r] 



(20b) 

On the quantum level, the evolution dictated by Hl()|) 
can be written as a map connecting the state of the sys- 
tem before each consecutive kick as 

|'0)„+l = UhUk\lp)n 

where Uh and Uk are, respectively, the evolution opera- 
tors for the harmonic oscillator and for the kicks. 



B. Open system: influence of environment 



d = mi>/2hx + i\Jl/2h'mvp^ 



(15) 



the above Hamiltonian can be written as 

OO 

H = hva)a + hKq cos [77(0 + a))] ^ 5{t - nr) , (16) 



where 



k\/h/2mv 



(17) 



is a scaling parameter related to the macroscopic limit - 
the so-called Lamb-Dicke parameter j^Oj. This parame- 
ter can be expressed as 77 = 27rAxo/A, where Axq is the 
width of the ground state of the harmonic oscillator, and 
A = 271 /k is the wavelength of the kicking force. This is 
the classicality parameter for the model under consider- 
ation. Its square is seen to be the ratio between h and 
the action 2mv/k'^ , playing the role of the dimensionless 
parameter fie// mentioned in the Introduction. The limit 
77 ^ can be achieved by letting fc ^ 0. This param- 
eter can be easily changed in ion trap experiments, by 



The influence of the environment on the system can be 
described by the master equation 



dp 
'dt 



H,p 



Cp, 



(22) 



where p is the reduced density operator of the system 
in the interaction picture. The first term on the right 
hand side of Eq. (|22|l corresponds to the unitary dynam- 
ics while the second term represents the non-unitary ef- 
fect of the environment in the Lindblad form 



Cp: 



{-fi/2) ( 2 Cj p cj ~c\cip~'p c\ Ci 



(23) 



The operators q are related to the form of system- 
environment couplings and the constants 7^ measure the 
strength of these couplings. 

Equation (|23|l is frequently found in the description of 
dissipative systems. It can be derived under very general 
assumptions, namely Markovicity and complete positiv- 
ity of the time-evolution of the reduced density operator 
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of system |31|, |33 • This last condition is defined in the 
following way. Let A be the system for which the re- 
duced density operator is defined, Ha the corresponding 
Hilbert space, and Aa the time-evolution map for the 
reduced density operator p. Consider any possible ex- 
tension of Ha to the tensor product Ha ^ Hb, where 
Hb any arbitrary Hilbert space; then is completely 
positive on Ha if (g) Ib is positive for all such exten- 
sions. Complete positivity corresponds to the statement 
that, if system A evolves and system B does not, any 
initial density matrix of the combined system evolves to 
another density matrix. 

In trapped ions one can use the techniques of "reser- 
voir engineering" [^|23,IS3 to build different kinds of Ci 
operators for the center-of-mass motion of the ion, and, 
in particular, the dissipative and diffusive reservoirs dis- 
cussed previously in the context of classical dynamics. 



1. The dissipative case 

Dissipation by a zero-temperature reservoir, in the 
rotating-wave approximation, is described by 122|l using 
just one operator ci = a 



dp 
'dt 



(2 



a pa' 



a) a p 



pa) a) 



(24) 



where H' has the same for as H given by (|16|) . with v 
replaced by the frequency given by Hll)|l . We will see 
that, with this choice, the oscillation frequencies of the 
quantum and classical systems will coincide. 

From the master equation H24I) , one gets the equations 
of motion for the expectation values between the kicks: 

{h)^Tr{af)^-i^l{d)-^{a), (25) 



(26) 
(27) 



-iVl{a) 

which can be written in terms of x and p as 
m 2 



One should note that, differently from the classical 
equations of motion, dissipation appears here in a sym- 
metric way with respect to position and momentum. This 
is related to the rotating-wave approximation, adopted in 
deriving H24|l : indeed, this approximation requires that 
the oscillator suffers many oscillation within the decay 
time (that is, one should have V <^v), which implies that 
the effect of the dissipation gets distributed between the 
canonical coordinates. 

Taking the derivative of H26|l and using H27II one gets: 



VL^ + — 



Tlx) = 0. 



(28) 



Using (|10l) . we can see that this equation, that describes 
the quantum dynamics between the kicks, is identical to 
its classical version so that the quantum and the 
classical systems oscillate with the same frequency. 



2. The diffusive case 

The purely diffusive reservoir master equation can ob- 
tained from H22|l by choosing two operators ci = a and 
C2 = with the same rate 71 = 72 = 7 



{2dpa) - "tAA_ ;,^t; 



a' ap — pa' a) 



+ (2a) pa — aa) p — pad)) 



(29) 



This is a combination of cooling and heating reservoirs 
and due to the fact that they have the same rates, all 
terms leading to drifts are canceled out and only diffusion 
terms survive. This fact becomes clear when one writes 
explicitly the corresponding Fokker-Planck equation for 
the Wigner function: 



dW 
~dt 

or, in terms of v and w. 
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dada* 



dW 
'df 



W. 



(30) 



(31) 



This equation is equivalent to the third term of 
rewritten in terms of the rescaled variables if we set 
7 = V/rf. 

A purely diffusive reservoir can be produced by random 
electric fields [s^ and is known to model the heating 
of vibrational energy observed in recent experiments on 
ion dynamics [23 |. 



IV. RESULTS 

The classical description of a chaotic dynamical sys- 
tem, either using single trajectories or a probability dis- 
tribution, is based on the analysis of the phase space and 
its structures. The definition of a single trajectory in 
the quantum case is prevented by the uncertainty prin- 
ciple and the suitable description of the system is based 
on quasi-probability distributions. The Wigner function 
fulfills almost all the requirements for being a true proba- 
bility distribution, as it is the only quantum distribution 
that yields the correct marginal distributions for any di- 
rection of integration in phase space, however it can ex- 
hibit negative values. For our purposes, this turns out 
to be an advantage, because it highlights the differences 
between quantum and classical dynamics. As a matter of 
fact, it is much easier to detect quantum signatures with 
the Wigner distribution rather than with the Husimi or 
Q-function. 

Oscillations between negative and positive values in 
the Wigner function are a sign of existence of quantum 
interference phenomena, absent in its classical counter- 
part. The role played by decoherence in washing out in- 
terference patterns is also easily visualized in the Wigner 
function \\^ . More than a visualization tool, the Wigner 
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function can be useful to derive some analytical results 
concerning the quantum-classical limit. 

In what follows we make use of the Wigner function 
and of its Fourier transform, the characteristic function, 
to obtain new results concerning the time scales for the 
quantum-classical transition. Combining the interaction 
with the environment and the possibility of varying the 
effective Planck constant, we are able to discuss not only 
the regions of parameters for the classical limit but also 
the behavior of the breaking time in open systems. 



A. System without reservoir 

In the absence of the interaction with the environment, 
the classical limit is investigated by changing the scaling 
parameter rj. One should remark that, in terms of the 
variables u and v, the initial distribution does not depend 
on ry, and is taken to be the same for the classical and 
the quantum systems. Of course, in terms of the original 
variable u and v, decreasing rj leads to a shrinking width 
of the initial distribution, both in the classical and the 
quantum situation. In any case, changing rj will affect 
both the classical and the quantum solutions, since the 
initial states are always taken to coincide, and a broader 
initial packet will explore, since the beginning, a larger 
region of phase space. 

In Fig.|Slwe show the Wigner function for 77 = 0.5 (top) 
and ?7 = 0.1 (bottom) after nine kicks, corresponding, re- 
spectively, to the classical situation depicted in Fig. Eli 
and its scaled version (not shown). Decreasing the value 
of the effective Planck constant one gets a quantum phase 
space that resembles more and more the overall classical 
structure but still with the presence of interference pat- 
terns. 

After some time these differences between quantum 
and classical evolutions become important and estimates 
can be made using the characteristic function C(A, A*), 
defined as 



C(A,A*) = Tr 



pe 



(32) 



Expanding the exponential of the cosine function in 
the quantum map (|21(l in terms of Bessel functions, it 
is possible to write explicitly the characteristic function 
after the n-th kick in terms of its initial value, Co(A„, A* ), 
as [13 

00 

C„(A, A*)= ^ Jrni{zi)Jm2{z2) ■ ■ ■ Jmr,iZn)Co{Xn,Ki)i 

(33) 



nil — — 00 



where 



Afe = Afe_ie'" + imfcTy, 
K 

Zk = 2Kq sin(^fe) = — sin(^fc) 
-^(Afe + A^), 



Ao 



A, 



(34a) 
(34b) 

(34c) 
(34d) 



u 0- 



u 




FIG. 5: (Color online) Wigner distribution after 9 kicks for 
•q — 0.5 (a) and rj = 0.1 (b). In both cases the Wigner function 
presents negative values but as the Lamb-Dicke parameter 
is decreased (closer to the classical limit) there is a better 
correspondence with the overall classical structure shown in 
Fig.Hd. 



J„, are Bessel functions and a ~ 



as in the classical 



It is interesting to compare this expression with the 
corresponding one for the classical system, which can be 
obtained by introducing the classical map (|20|l into the 
appropriate classical definition of the characteristic func- 
tion. This definition follows from the quantum expres- 
sion H32|) upon changing the trace by a double integral, 
the operators by complex numbers and the density ma- 
trix by a classical probability density. One gets then: 

1 



J2 -^'"1 (5^1 

mi ,. . . ,m„ — — CX2 ^ 



C(^i)(A,A*) = 

mi 

X Co(A„,A*) 



K 

■ Jrrin 1 ~9 



(35) 



It is clear that the classical expression is obtained from 
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the quantum one when \^k \ ^ 1, in which case 
sin(^fe) « ■ 



(36) 



Since, according to Eq. (|34a|) . is proportional to rj, 
this approximation should hold for sufficiently small ij in 
the beginning of the evolution of the system. However, 
as time evolves, and |Afe| grows, it eventually ceases to 
be true. This is precisely where the breakdown between 
classical and quantum evolution occurs. 

An expression for the breaking time can be obtained 
by comparing the quantum and classical characteristic 
functions. We may define it as the time at which the 
approximation H3fci|) fails, or, in other words, for which 

~ 1. Assuming the strong chaos condition {K ^1), 
one is able to derive [l^ 



Th 



\n{2K/r]) 
ln{K) ' 



(37) 



where K = Ksm{a). This result displays the scaling of 
the breaking time with the logarithmic of l/yy^, which 
stands for the tigj-f already mentioned in the Introduc- 
tion. A numerical check of this expression needs an oper- 
ational definition for the breaking time, which involves, 
also, a choice of an appropriate measure of the distance 
between quantum and classical systems. Information 
measures that can be used to compare two different dis- 
tributions are available in the literature [3l|, and have 
been applied in the context of quantum-classical tran- 
sition for chaotic systems jlj]. Measures based on the 
comparison between the whole distributions, although 
more complete, would lead to a large increase in com- 
putation time and to experimental difficulties. Although 
the Wigner function has already been measured in ex- 
periments with trapped ions [3.6[ , it would be challenging 
to resolve the details of the interference fringes seen, for 
example, in Fig. [S}b. The relative distance between the 
classical ((Aw^;)) and quantum ((Aw^)) variances of the 
distributions, defined as 



dr — 



(38) 
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FIG. 6: Classical and quantum variances (left) and the rela- 
tive distance d,- (right) as a function of the number of kicks 
for rj — 0.5 (top) and rj = 0.1 (bottom). The quantum (solid) 
variance remains close to the classical (dashed) for a longer 
time when the Lamb-Dicke parameter is smaller. The break- 
ing time, indicated by the arrows, correspond to the instant 
at which the relative distance gets larger than a chosen value 
e — 0.1 (horizontal line in the right panels). 
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FIG. 7: Breaking time as a function of ln(l/r;) for K — 2.0 
and 5 = 6. Despite the oscillations there is a clear linear 
behavior confirming the scaling predicted by (^JJ. 



is a much simpler quantity that already shows the scaling 
properties of H37|l . Fig.|H|shows the classical and quantum 
variances (left) and the relative distance d,- (right) as a 
function of the number of kicks for two different Lamb- 
Dicke parameters. The separation time is defined as the 
time at which the relative distance crosses a given value 
e (e = 0.1 in the figure). In Fig. [7|we plot r^j obtained 
in this way as a function of ln(l/77) and, although the 
absolute value of the breaking time depends on the choice 
of e, tests with e ranging from 5% to 30% show only 
slight modifications in the curves and confirm the scaling 
behavior (|37|l . independently of the particular definition 
of the separation. 



B. Dissipative environment 

The analytical solution for the quantum dissipative 
problem in terms of the characteristic function is also 
given by H33|) and the only change affects the rela- 
tion (|34a|l that must be replaced by 

Afc = Afc_ie'"e-^"/2 + imkT], (39) 
1 

where F was introduced in the classical case. Besides the 
usual rotation due to the harmonic motion represented by 
the complex exponential in (|39|) there is also an exponen- 
tial decay due to the dissipative drift in the characteristic 
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TABLE L Breaking time for dissipative dynamics in the four 
different regions of parameters rr/2, 77 and K' . 



2.0 



Nonlinearity 


Deep quantum 


Weak quantum 


1.5 


strengtfi 


regime 


regime 






rr/2 < ln(r;/2) 


rr/2 > ln(ry/2) 


dis 

^1.0 


in (7?') > rr/2 


(a) rh ~ 1 kick 


(„\ ^ ^ ln(2if'/,,) 
(^1 'h ~ in(K')-rT/2 


0.5 


In (7?') < rr/2 


(b) rfi 1 kick 


(d) Th ^ 00 





function's argument |37| . 

In analogy with what was done for the system without 
reservoir, one can obtain the breaking time by examining 
when the quantum characteristic function cannot be de- 
scribed anymore by its semiclassical approximation. This 
procedure is fully described in Appendix El and the re- 
sults are summarized in Tabled The first column repre- 
sents a region characterized by a deep quantum regime 
where a classical description of the system is already not 
valid right after the first kick, for all finite values of the 
dissipation. It is interesting to note that, in this region, 
quantum-classical correspondence is lost even if dissipa- 
tion is sufficient to bring the classical system into a peri- 
odic regime. 

In the second column we have the most interesting 
range of parameters concerning the quantum-classical 
transition (we call it "weak quantum regime" ) , where 
two different regimes exist: one indicating an increase of 
the breaking time with dissipation, and another showing 
close quantum-classical behavior for all times {Tn — s- 00). 
This latter case corresponds to a situation where the dis- 
sipation is so strong that classical chaos is suppressed 
and the system goes to a simple attractor. The breaking 
time for region (c) in Tabled 



dis 



In {2K'/r]) 
In (i?') - rr/2' 



(40) 



with K' — K'sm(a), increases as the dissipation coeffi- 
cient grows, but it keeps the same logarithmic-scale de- 
pendence with respect to the effective Planck constant as 
in the case without reservoir. 

Although r^*^ could be arbitrarily large, as pointed 
out by lomin and Zaslavsky 'S| in a recent derivation 
of an expression similar to H4Q(I . this is not the case if 
one wants to preserve a strange attractor. The condition 
rr/2 = \n{K'), which separates regions (c) and (d) in 
Table U corresponds to the situation where the origin of 
phase space changes from an unstable to a stable fixed 
point. However, instability of the origin is not sufficient 
to ensure chaotic dynamics. Indeed, Fig. |31 exhibits a 
large range of values of rr/2 for which the system is 
attracted to some periodic trajectory, even when rr/2 < 
In {K). 

The ratio between the breaking times with and without 
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riG. 8: Ratio between tfie breaking time witfi and without 
dissipation as a function of rr, for the same parameters as 
in rig. 1^ and 77 — 0.5. Solid line shows the analytical predic- 
tion I41I I while triangles show numerical simulations compar- 
ing quantum and classical systems. The differences between 
quantum and classical environments are responsible for the 
deviations observed for larger values of rr. Smaller values of 
rr where not considered due to numerical constraints. 



dissipation is given by 



Th 



In [K') 



In (if') - rr/2 



(41) 



This expression exhibits the increase of the breaking time 
as a function of rr. In Fig. |S1 this relation is plotted to- 
gether with numerical simulations for K' = 6.0 with the 
horizontal axis ending at the value (Tt /'2,)max ~ 0.51 of 
the dissipation constant, for which the Lyapunov expo- 
nent becomes negative (see Fig.O. For this value of the 
kicking strength the maximum increase in the breaking 
time is around 1.5 and therefore does not help consider- 
ably to achieve the classical limit. This increase depends 
on the values of the chaoticity parameter but even for 
very large values [K' k, 500) it is less than a factor of 4. 

Although expression (|41|l is independent of the scal- 
ing parameter 77, some remarks about the role played by 
the macroscopic limit in the dissipative case are needed. 
It should be noted that the quantum description of a 
zero-temperature reservoir given by (|24|l is not com- 
pletely equivalent to the classical description based on 
the map H12|l . In fact, a classical distribution submitted 
only to the dissipative dynamics would contract towards 
the origin, while a quantum distribution would end up 
in the ground state, which has a finite width due to the 
uncertainty principle. This argument does not invalidate 
the results presented in Table ^ but rather emphasize 
that they are valid for the semiclassical approximation 
that is not equivalent to the fully cl^issical system based 
on map (|12|1 . 

This discussion also provides an explanation for the de- 
viations between numerical and analytical results shown 
in Fig.|Sl Only in a small range of the dissipation strength 
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(0.12 < Vt/2 < 0.22) we were able to see the expected 
growth in the breaking time, while for larger values of 
Ft/ 2, the effects of the difference between classical and 
quantum systems, discussed above, become predominant. 

It is interesting to note the sudden growth in the break- 
ing time for the region corresponding to the periodic win- 
dow around Tt/2 = 0.4 shown in Fig. El This should be 
expected, since the quantum and classical systems should 
stay together for a longer time outside the chaotic region. 
On the other hand, as Fr increases, inside the same re- 
gion, one notices that the breaking time decreases. This 
is due to the fact that, for larger dissipation, the distribu- 
tion shrinks at a faster rate, implying that the two distri- 
butions approach at an earlier time the region around the 
origin, where the uncertainty principle plays a dominant 
role. 

For smaller values of the dissipation parameter, the 
system could spread over a large region of phase space 
demanding a huge amount of computational resources. 
This, together with reliability problems for even smaller 
values of Ft/2, imposed the limit Ft/2 > 0.12 for the 
dissipation parameter in our calculations. 

It is important to understand the meaning of the 
breaking time and its consequences for the dynamics of 
the system at different times. In particular, the station- 
ary state produced by the dissipation is of much interest 
and this issue of long time behavior has been addressed 
before in the case of the standard map |^ . The existence 
of a finite Tri means that quantum and classical dynamics 
cannot be equivalent for all times but does not necessar- 
ily mean that they have to be different for all t > th. In 
fact, the numerical simulations for the evolution of the 
variance of the dissipative KHO show that, in some cases, 
quantum and classical calculations share the same final 
stationary behavior but with different transient regimes 
as can be seen in Fig. |5| One should note, however, that 
this is not necessarily true for the whole phase space dis- 
tributions, which can be different, although having the 
same second moments. This can be illustrated through 
the comparison of the Wigner function depicted in Fig. 1101 
and the strange attractor shown in Fig.^ The quantum 
distribution clearly does not exhibit all the structures 
presented in the classical case despite the fact that it lies 
in the region that contains the classical attractor. 



C. Diffusive environment 

The evolution of the system under the influence of 
the diffusive environment can be solved analytically in 
terms of the characteristic function, as described in Ap- 
pendix^ The solution between consecutive kicks can be 
written as 

C{X,X*,t) = C(A,A*,0)e-^l^l'^. (42) 

Using this solution together with , we can establish 
a recurrence relation for the characteristic function after 
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FIG. 9: Classical (dashed line) and quantum (solid line) evo- 
lution of (Aw^) as a function of the number of kicks for 
K' = 6.0 and rT/2 = 0.36. Classical and quantum variances 
share the same asymptotic behavior despite the fact that they 
can differ in a transient regime. 
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FIG. 10: (Color online) Wigner function for the same param- 
eters of the strange attractor of Fig.0|and rj = 0.5. The quan- 
tum distribution, although lying in the region of the classical 
attractor, does not show the classical small-scale structures. 

the n-th kick as 

oc 

C„(A, A*) = 6-^1^1'^ J2 Jmi(^i)C„-i(Ai,At), (43) 

m 1 — — oo 

with all the variables defined as in (|34|) . 

While in the dissipative case the drift effect adds an ex- 
ponential factor to the arguments of the Bessel functions, 
in this case, a Gaussian multiplies the whole sum. This 
difference is crucial for understanding the influence of 
this environment in restoring the quantum-classical cor- 
respondence. First, one should note that the role played 
by the Lamb-Dicke parameter in the above expressions is 
the same as in the case of a system without reservoir and. 
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therefore, the macroscopic approximation would lead to 
the same result as before. It is clear, however, that dif- 
fusion should have an important effect on the behavior 
of the system. This can be seen through a more careful 
analysis of expression H43() . 

Assume that the quantum and classical dynamics co- 
incide at kick n and forget, for the moment, the diffusion. 
As discussed before, the two dynamics will differ as long 
as the approximation sin(^) « ^ fails. This gives an esti- 
mate of the values of that lead to quantum corrections 
as 



Id 



i(A. 



> 1. 



(44) 



This equation shows that quantum corrections are asso- 
ciated with large values of A and we can define the typical 
values At for which the corrections appear as 



> 



V 



(45) 



The effect of diffusion is to cut off the contributions 
from large values of A, due to the Gaussian modulation 
in 142|) . This implies that the values of A that satisfy 145|) 
may be attenuated by the Gaussian pre-factor, which ren- 
ders them inefficient in promoting the quantum-classical 
separation. 

This may become more intuitive if we go back from the 
characteristic to the Wigner function: because they are 
related by a Fourier transform, the larger values of the 
characteristic function correspond to the small-scale in- 
terference structures in the Wigner function and thus to 
quantum corrections. The disappearance of these small- 
scale structures, due to the Gaussian modulation of the 
characteristic function, has been extensively discussed in 
the literature [H IHIll El 13: it leads to a better 
correspondence between quantum and classical distribu- 
tions, and to the emergence of the classical world [38l |. 
Fig. II II shows, for two different Lamb-Dicke parameters, 
classical and Wigner distributions in the presence of dif- 
fusion. These distributions are much more similar to each 
other than the corresponding distributions for the system 
without reservoir, displayed in Fig. [S] The importance 
of diffusive effects will depend, however, on the value of 
the effective Planck constant, as can be seen by com- 
paring the impressive correspondence between quantum 
and classical distributions for ij = 0.1 (right) and some 
evident differences that persist for rj = 0.5 (left). 

Although it is hard to define precisely the value of the 
diffusion needed to restore the classical limit, a rough 
estimate can be obtained as follows: if the values of At, 
given by (|45|l lie outside the range defined by the width of 
the Gaussian, then the corrections should remain small. 
Using H45|l . this conditions reads 



(46) 



where D = 7t/2 plays the role of a dimensionless dif- 
fusion coefficient for the renormalized coordinates {u,v). 




FIG. 11: (Color online) Classical (top) and Wigner (bottom) 
distributions for K = 2.0, D = 0.01047 and ?? = 0.5 (left) or 
77 = 0.1 (right). Diffusion leads to a better quantum-classical 
correspondence as compared to Fig.|3 For — 0.1 this corre- 
spondence is quite impressive while for ri — 0.5 differences still 
remain. Diffusion also prevents the appearance of small-scale 
structures on classical dynamics (top). 



This simple argument shows that there is a critical dif- 
fusion for which the classical and quantum dynamics re- 
main close to each other, and that it scales as 



Dci- oc 77 . 



(47) 



In terms of the diffusion coefficient in Eq. H13|l . one 
has: 



Per OC Ti'^/t . 



(48) 



This is the diffusion coefficient that corresponds to the 
non-renormalized variables v and u. This result is consis- 
tent with those found for example in Refs. [U and 
One should expect, however, that the strength of the non- 
linearity, represented by K in our case, should play an 
important role in such a scaling law. The argument lead- 
ing to H47|l was based in the estimation of the values of ^1 
at which quantum corrections become important, with- 
out taking into account the size of these corrections. We 
have not studied in detail the actual separation between 
the two distributions, as they evolve with time. This is 
the reason why our simple argument could not account 
for the influence of nonlinearity, which is hidden in l|47f) . 
A detailed investigation of the separation time for the dif- 
fusive case, or an estimation of thel error introduced by 
the neglected contributions, would certainly display this 
dependence. Scaling relations between effective Planck 
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FIG. 12: Relative difFerence between quantum and classical 
variances as a function of the number of kicks for = 0.5 
and D equal to (solid); 0.00209 (dashed) and 0.0209 (long- 
dashed). For large enough D, the relative difference remains 
always smaller than the percentual threshold e represented by 
the horizontal lines. 



constants, environment and nonlinearity strengths in the 
context of the quantum-classical transition have been ob- 
tained by many authors 0)0,O)E)0| and have moti- 
vated recent interest [l4'| in finding the properties of such 
scaling. 

The above considerations suggest that the breaking 
time should diverge when the diffusion coefficient exceeds 
a certain critical value Dcr- This may be easily under- 
stood from Fig. ^1 which displays the time evolution of 
the relative distance for the quantum and classical vari- 
ances, for several values of D. One should note that, as 
D becomes larger than a critical value, which depends on 
the percentual threshold e adopted for the definition of 
the separation time, the relative distance remains always 
smaller than this threshold, implying an infinite separa- 
tion time. On the other hand, for sufficiently small diffu- 
sion coefficients, one should recover the logarithmic time 
scale. Although we were not able to derive an analyti- 
cal expression for the breaking time when D < Dcr, our 
numerical simulations show that it remains indeed prac- 
tically identical to the result obtained when no reservoir 
is present, in this region of parameters. 

The behavior of the breaking time in presence of a dif- 
fusive environment is shown in Fig.^J For D = 0.01047 
(top), one can see that the breaking time basically lies on 
the curve corresponding to the system without reservoir, 
and increases abruptly at a given value of the Lamb- 
Dicke parameter (rj = 0.31, indicated by an arrow in the 
figure), with some of the oscillations, already present in 
the D — case, amplified. For rj < 0.31, the breaking 
time experiences a sudden growth 4 the corresponding 
points are not shown in the figure (numerical tests were 
performed for the maximum of 50 kicks up to 77 = 0.2). 
For r] = 0.5 (bottom), the increase in the noise introduces 
very small changes in the breaking time, which grows 
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FIG. 13: Top: breaking time as function of ln(l/r7) for D — 
0.01047 (triangles) and D — (line). The two curves are 
essentially the same until rj ~ 0.31 (shown by the arrow) 
when there is no separation anymore. Some peaks of the 
D — curve are amplified, due, probably, to the definition of 
breaking time adopted. Bottom: breaking time as function 
of D for rj = 0.5. The separation time remains basically 
the same, with a small increase (« 10% at maximum), and 
suddenly increases for Dc- ~ 0.13. 



quickly when D w 0.01466. Again, the abrupt increase 
indicates that the differences between quantum and clas- 
sical variances remain bounded below a given limit e. 

One should note, however, that the small changes in 
the breaking time for D < Dcr do not imply that the 
environment has no effect at all in the dynamics. In fact, 
observing again Fig. 1121 we see that the maximum dis- 
tance between the variances decreases smoothly with in- 
creasing noise strength and may eventually come to zero, 
indicating a perfect quantum-classical correspondence. It 
is interesting to note how these twoj^different quantities, 
breaking time and maximum distance, give complemen- 
tary information about the dynamics. 

From these results one could infer, naively, that diffu- 
sion is sufficient to restore the classical limit for a chaotic 
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system no matter the value of the effective Planck con- 
stant. Indeed, for a given Lamb-Dicke parameter, one 
can always find a large enough diffusion coefficient to 
bring the quantum and classical dynamics sufhciently 
close to each other. Nevertheless, such a statement de- 
serves some reserve. In fact, diffusion washes out not only 
the interference pattern in the Wigner function but also 
the structures in the classical phase space (see Fig. Illll . 
and one could claim, therefore, that when the critical dif- 
fusion coefficient is very large, the chaotic characteristic 
of the system is lost and the system follows, basically, a 
diffusive dynamics induced by the environment. A simi- 
lar situation was described in the dissipative case where a 
large enough dissipation was sufficient to suppress chaos 
and bring the system to a periodic behavior. There, how- 
ever, we could clearly distinguish between chaotic and pe- 
riodic behavior through the calculation of the Lyapunov 
exponent, while here, even though generalizations of Lya- 
punov exponents for distributions exist 39, 40], there 
is no sharp distinction between chaotic and regular be- 
havior. The description of the system can become even 
more complicated with the addition of diffusion in view of 
the mixed phase-space structure of the system. Indeed, 
the difficulties in characterizing chaos lie not only on the 
smoothing of the phase space structures (see Fig. ITTla'l. 
but also on the fact that the distribution can flow from 
regular to chaotic regions that were well separated when 
no reservoirs were taken into account. 

Finally, one should remark that we have not analyzed 
here the long-time behavior of the system. In particular, 
it is well-known that diffusion may affect quantum local- 
ization, which occurs at times much longer than those 
here considered 0, H, IS Q • In the absence of a reser- 
voir, the long-time dynamics of the model analyzed in 
this paper may display either a quantum diffusion or a 
baUistic behavior Ha. 



V. CONCLUSIONS 

We have shown that it is possible to discuss separately, 
in a physically relevant way, the roles of the macroscopic 
limit and of different system-environment interactions in 
the quantum-classical transition for a chaotic system. We 
have considered the kicked harmonic oscillator coupled 
to two distinct reservoirs, giving rise in the classical limit 
to either pure dissipation (zero-temperature reservoir) or 
pure diffusion (random force), in a situation that could 
be implemented in state-of-the-art ion trap experiments. 

In the chaotic regime, when the interaction with a 
reservoir is not taken into account, the classical and quan- 
tum dynamics start diverging after a time that depends 
logarithmically on the ratio between a typical action of 
the system and Planck constant. We have used an opera- 
tional definition of the breaking time in terms of measur- 
able quantities, which would allow an experimental test 
of this logarithmic time scale. 

In the dissipative case we established regions of pa- 



rameters corresponding to different time scales. There 
is a region where quantum corrections appear right af- 
ter the first kick and quantum-classical correspondence 
is lost already in the beginning of the evolution. De- 
creasing the Lamb-Dicke parameter, one reaches a re- 
gion where quantum-classical correspondence persists for 
a time that, as in the system without reservoir, grows 
only logarithmically with the classicality parameter. We 
have also shown that, for a fixed effective Planck con- 
stant, close agreement between quantum and classical 
predictions is only possible for dissipation strengths large 
enough to bring the system into regular behavior. 

In the diffusive case, we were able to establish that the 
behavior of the quantum-classical separation should be 
markedly different, depending on whether the diffusion 
coefficient is above or below a certain critical value Dcr- 
For D > Dcr, this separation should remain small, and 
infinite separation times may even be obtained, at val- 
ues of D that depend on the definition adopted for the 
critical percentual separation. We have also presented 
numerical evidence that, for diffusion coefficients below 
this limit, the breaking time behaves as in the case with- 
out reservoir. Furthermore, we obtained an analytical 
estimation of the dependence of the critical diffusion co- 
efficient on the effective Planck constant, which shows 
that the farther away from the classical limit is the sys- 
tem, the larger must be the effect of the environment to 
restore the quantum-classical correspondence. 

Although the coupling with the environment helps re- 
store the quantum-classical correspondence for a system 
that is close to the macroscopic regime, for systems in a 
deep quantum region the critical diffusion can be so large 
that it brings classicality at the expense of reducing, or 
even extinguishing, the chaotic features of the system. 

The behavior of the system under the influence of other 
kinds of environment could also be explored in this con- 
text. Thermal and phase reservoirs are examples of dif- 
ferent environments, already produced in ion trap exper- 
iments, that could be used for this purpose, and could 
lead to interesting results for the quantum-classical tran- 
sition scenario. 
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APPENDIX A: BREAKING TIME 

It is not difficult to generalize the solutions 
and H34|) , which were obtained for the characteristic func- 
tion corresponding to a system without reservoir, to 
the dissipative regime. Between two consecutive kicks 
the harmonic evolution just rotates the system in phase 
space, and this effect appears only in the complex expo- 
nential in H34a|) . The solution of the dissipative master 
equation simply adds an exponential decay e"'"'^/^ 
and, therefore, the solution of the full problem is still 
given by (p?!^ : 

oc 

^^(A, A*)=: Jmi (21) Jm2 (^2) • • ■ '/jTln (2:71) C'o(A„, A* ), 

mi ,...,m,i — — 00 

(Al) 

where 

Afe Afe_ie*"e-^^/^ +i2mfc77, (A2a) 
K' 

Zk = 2Kq sin(a) = sin(a), (A2b) 

a = -f(Afc + A:), (A2c) 
Ao = A. (A2d) 

The procedure to obtain the breaking time in the dis- 
sipative case follows very closely the one used in 18| for 
the situation without reservoir. First we should note that 
the macroscopic limit, as discussed previously, is achieved 
by letting ry — > 0, which means that the sine functions 
in ljA2b|l can be approximated by their argument, i.e., 

sinfe)^^. (A3) 
In this limit, we obtain the scmiclassical characteristic 



function 

00 

C„(A,A*) = J,nA2K,^l)...JmA'2Kq^n) 

mi ,.. .,mn — — 00 

X Co(A„,A:), (A4) 
with initial condition 

/oo 
dV^'o(/i,Ai*)e^"''*-'>. (A5) 
-00 

When no reservoir is present, the semiclassical approxi- 
mation leads to the classical characteristic function de- 
rived, directly, through the corresponding classical map. 
Here this is not the case, due to the fact that the zero- 
temperature reservoir leads to distinct features in the 
quantum and the classical models. If submitted only to 
the dissipative dynamics, the quantum system will end 
up in the its ground state, which has a finite width, while 
a classical probability distribution would contract to a 
point located at the origin. In this way, the semiclassical 
characteristic function shows exactly the same nonlinear 
dynamics as the classical but with an intrinsic quantum 
property due to the uncertainty principle. Of course, this 
effect becomes smaller as quantum fluctuations become 
negligible compared to the size of the system, which oc- 
curs for small values of the Lamb-Dicke parameter. 

The replacement of (ESI by l|A4|l is only valid if 
holds for every k. Taking into account that the Bessel 
functions decrease exponentially for \mk\ ^ 2Kq^k, we 
can truncate the sums in (jAlfl . by estimating the maxi- 
mum values of jm^l in each sum. The relevant contribu- 
tions are the ones in which the Bessel function index is 
of the order of its argument and, therefore. 



\m2\ 



2t] 



-Trjl 



f\ i2a I \* -i2a\ -2rr/2 o™ „• f -\ -Tt/2 

(Ae + A e je ' — imirj sm(aje ' 
(Ae™" + A*e-"")e-"^^/2 _ 3^^^ sin[(n - l)a]e-("-i)^^/2 - ... - 2m„_i77 sin(a)e 

I 



-rr/2 



Considering, now, the strong chaos limit {K' ^ 1), we 
have 



K' 



mi 



2rj sin(a) 



e-^^/^...,|m„| 



K' 



2r] sin(a) 



-nrT/2 



ratio between two consecutive ^'s is given by 
while the first argument is 



From the above considerations we can study the regime 
of validity of the approximation <C 1 for all k. The 



la 



'Vrl2 



(A6) 



(A7) 



From ljA7|l we can get two different conditions on the 
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first argument |^i|. If Tt/2 < ln(?7/2) then |^i| > 1 
and, right after the first kick, the classical approxima- 
tion is not valid anymore, so quantum predictions should 
differ from the classical ones. This corresponds to the 
results shown in the first column of Table (QJ, which are 
represented by r^" « 1 kick. If, on the other hand, 
Tt/2 > ln(?7/2) and, therefore, |^i| < 1, then we have 
two new possibilities depending on condition ljA6|l . 

Choosing Tt/2 > ln{K'), the ratio between two con- 
secutive ^'s is less than one and the sequence of ^^'s is 
decreasing with increasing k. Because |fi| < 1, all the 
terms will be smaller than one and quantum and classi- 
cal evolutions should stay close to each other at all times 
(r^*'' oo). Nevertheless, if we have Tt/2 < \n{K') the 
sequence of ^^'s is increasing and one should expect that 
there is a for which the condition < 1 is not ful- 
filled anymore and a breaking time will exist. This will 
happen when w 1, which can be expressed as 



The diffusive dynamics is now introduced by the re- 
placement of p in the above equation by H29|l , which gives 



C(A,A*) = ^Tr 



-a^ape^^'-^'^ + 2apa) e^^' '^'^ 



+2a''pae^^'-^'^ - paa^ e^^' - ^' ^] . (B2) 



Rewriting the exponentials using the Baker-Hausdorff 
formula and the ordering properties 



e-''°V(a,a^)e'^"' = /(a + /3,o^), (B3a) 
e''V(a,a^)e-^' = /(a,aU/?), (B3b) 



-nTT/2 



2K' 



1. 



(A8) 



we obtain 



Taking the logarithm and noting that n, the number of 
kicks, corresponds to the time in units of r, we arrive at 
an estimate for the breaking time as 



\TL{2K'/ri) 
In(X') - Tt/2' 



(A9) 



C(A,A* 



-Tr 



~2|Ape^°'-^'V 



-r|ApC(A,A*). 



(B4) 



APPENDIX B: DIFFUSIVE CASE 



This equation can be readily integrated, giving as so- 
lution 



The evolution of the symmetric-ordered characteristic 
function defined in (|32|l is given by 



C(A,A*,i) = C(A,A*,0)e~^l^l'*. 



(B5) 



C(A,A*) = Tr pe 



^\a^ —A* a 



(Bl) 
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